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^ , ABSTRACT 

■ We describe new Chandra spectroscopy data of the cluster which harbors the proto- 
typical "head tail" radio galaxy 3C 129 and the weaker radio galaxy 3C 129.1. We 

{SJ . combined the Chandra data with Very Large Array (VLA) radio data taken at 0.33, 

' 5, and 8 GHz (archival data) and 1.4 GHz (new data). We also obtained new HI ob- 

C*"** I servations at the Dominion Radio Astrophysical Observatory (DRAO) to measure the 

CNI ■ neutral Hydrogen column density in the direction of the cluster with arcminute an- 

gular resolution. The Chandra observation reveals extended X-ray emission from the 
radio galaxy 3C 129.1 with a total luminosity of 1.5xl0 41 erg s _1 . The X-ray excess 
is resolved into an extended central source of ~2"(1 kpc) diameter and several point 
sources with individual luminosities up to 2.1x10 erg s _1 . There is no evidence for 
a correlation between the 3C 129.1 X-ray and radio morphology. In the case of the 

■ radio galaxy 3C 129, the Chandra observation shows, in addition to core and jet X-ray 

emission reported in an earlier paper, some evidence for extended, diffuse X-ray emis- 

q • sion from a region east of the radio core. The 12"x36"(6 kpc x 17 kpc) region lies "in 

!— i ' front" of the radio core, in the same direction into which the radio galaxy is moving. 

\ We use the radio and X-ray data to study in detail the pressure balance between the 

■ non-thermal radio plasma and the thermal Intra Cluster Medium (ICM) along the tail 

of 3C 129 which extends over 15'(427 kpc). Depending on the assumed lower energy 

• i-h , cutoff of the electron energy spectrum, the minimum pressure of the radio plasma lies 

1 a factor of between 10 and 40 below the ICM pressure for a large part of the tail. We 

discuss several possibilities to explain the apparent pressure mismatch. 

CS ' 

1 INTRODUCTION studied the radio plasma/ICM pressure balance for FR I and 

FR II radio galaxies, in most cases finding the minimum 

,. „ t, „ , .... , synchrotron pressure lower than the ICM pressure by ap- 

(Jbservations ot Radio Calaxies within Calaxy ■' r r j r 

Clusters ' proximately one order of magnitude (e.g. Feretti et al. 1990, 

Hardcastle & Worrall 2000 and references therein) . Further- 
Tailed radio galaxies within galaxy clusters have drawn in- more, the long tails of radio galaxies carry the imprint of 
tensive study for several decades on both theoretical and ob- relative motion between the non-thermal plasma and the 
servational grounds. Radio and X-ray data complement each ambient hot gas. In the parlance of the field, they reflect the 
other in that the radio data show high resolution images of weather conditions in the ICM. The potential of such obser- 
the jets and more extended synchrotron features yielding in- vations to reveal details of cluster mergers such as subsonic 
formation on strength and orientation of the magnetic field and transonic bulk flows, shocks, and turbulence has been 
within the radio plasma and between the radio galaxy and investigated on theoretical grounds by a number of authors 
the observer, while the X-ray data gives information about (Loken et al. 1995, Burns et al. 2002). On the observational 
the temperature, chemical composition, density, and pres- side, Gomez at al. (1997) studied a large number of Wide 
sure of the Intracluster Medium (ICM). Thus, the jets prop- Angle Tail (WAT) radio galaxies inside Abell clusters. They 
agate into a well defined medium, which greatly facilitates found a strong correlation between the orientation of the ra- 
the interpretation of the radio images. It becomes possi- dio tails and the elongations of the ICM surface brightness 
ble to study, for example, the pressure balance between the profiles, and interpreted it as evidence for bulk plasma mo- 
radio-plasma and ICM, and thus to search for hidden pres- tion that follows cluster mergers. Bending of radio tails has 
sure components in the radio plasma. Various authors have often been explained by invoking motion of the radio galaxy 
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in combination with "stormy ICM weather" (e.g. Feretti et 
al. 1985, Sakclliou, Merrifield and McHardy 1996) but this 
interpretation does not work for all bent radio tails (e.g. 
Eilek et al. 1984). 

In addition to facilitating studies of radio galaxies and 
the ICM themselves, observations of radio galaxies in galaxy 
clusters make it possible to study the interaction of the radio 
galaxies with the ICM and vice versa. The radio plasma can 
interact with the ICM and form so-called cavities devoid of 
ICM that show up as depressions in the X-ray surface bright- 
ness maps. Cavities have been observed in the Perseus clus- 
ter Bohringer et al. (1993), in the Cygnus A cluster by Car- 
illi et al. (1994); more recently, CHANDRA observations re- 
vealed cavities in the clusters Hydra A, A 2597, A 2052, and 
A 4059 (McNamara et al. 2000, 2001, Blanton et al. 2001b, 
Heinz et al. 2002). It is not yet clear whether these cavities 
are merely curiosities or, whether they influence the ICM 
in an important way by supplying a dynamically important 
amount of heat and/or magnetic field to the ICM (Briiggen 
et al. 2002, Quilis, Bower, & Balogh 2001, Reynolds, Heinz, 
& Begelman 2002, Fabian et al. 2002). The relative motion 
of the ICM and the Interstellar Medium (ISM) of the ra- 
dio galaxy's host galaxy can result in a modification of the 
ISM and in a bow shock forming "in front" of the radio 
galaxy (e.g. Stevens, Acreman, Ponman 1999, Toniazzo & 
Schindlcr 2001, and references therein). Although these in- 
teractions can also be observed in normal galaxies, the radio 
maps of radio galaxies can be used to better interpret the 
data by constraining the motion of the galaxies in the plane 
of the sky as well as the ICM magnetic field. 

The Cluster 3C 129 and Outline of the Paper 

In this paper we present Chandra data on the nearby galaxy 
cluster 3C 129 (z = 0.0223). Owing to its low Galactic lati- 
tude (I = 160° .5, b — 0.3°) it has not been studied extensively 
at optical wavelengths. In addition to velocity measurements 
of 3C 129 and 3C 129.1 (Spinrad, 1975), we are aware of only 
one other velocity-verified cluster galaxy: WEIN 048 (Nils- 
son et al. 2000). Optically, the galaxy 3C 129 has been de- 
scribed as a weak elliptical galaxy without peculiar features 
(Colina & Perez- Fournon 1990). X-ray observations with the 
Einstein, EXOSAT, and ROSAT satellites (Edge & Stewart 
1991; Leahy & Yin 2000) yielded a cluster luminosity of 
2.7xl0 44 ergs s _1 in the 0.2-10 keV energy band, a central 
gas density of 6xl0~ 3 cm~ 3 and an ICM mass of 3.6 x 10 13 
M©. Based on Einstein IPC data Leahy & Yin (2000) found 
evidence for a cooling flow with M « 85 M yr" 1 . However, 
reanalyzing the same data, Taylor et al. (2001) derived an 
upper limit of M < 1.2 M©yr~ . Based on the Chandra 
imaging spectroscopy data we performed a spectroscopic de- 
projection analysis, taking into account the asymmetry of 
3C 129's X-ray surface brightness by using a set of elliptical 
ICM shells (Krawczynski 2002, called Paper I in the follow- 
ing) . This analysis did not show a decrease of the ICM tem- 
perature toward the cluster center, as expected for a cooling 
flow. 

The cluster contains two radio galaxies. The radio 
source 3C 129 is a well studied galaxy in the radio wave- 
length range as a prototype head-tailed (HT) galaxy with a 
strongly curved, 15' long (426 kpc), two armed tail (Miley 
1973; van Breugel & Jagers 1982; Kassim et al. 1993; Fer- 



etti et al. 1998). In an earlier paper (Harris, Krawczynski & 
Taylor, 2002, called Paper II in the following) we presented 
evidence for X-ray emission from the inner northern jet of 
3C 129. Near the projected center of the X-ray galaxy clus- 
ter is the weaker and smaller radio source 3C 129.1. This 
galaxy also has dual radio jets which extend over ~2 arcmin 
(57 kpc) (Downes 1984; Jagers 1987; Kassim et al. 1993). 

In this paper we describe X-ray and radio observations 
of this cluster. We complement the Chandra X-ray data with 
archival 0.33 GHz, 5 GHz and 8 GHz and with new 1.4 GHz 
Very Large Array (VLA)* radio data. We use the radio and 
X-ray data to search for evidence of interactions between the 
radio plasma and the ICM and to study the pressure balance 
between the two plasmas. Since the cluster lies at low Galac- 
tic latitude we expect substantial variations of the neutral 
Hydrogen column density N H over the field of view of the 
Chandra observations. We therefore obtained HI observa- 
tions with the DRAO telescope. Tables 1 and 2 summarize 
the data sets used in our analysis. 

The remainder of the paper is structured as follows. We 
introduce observations and data sets in Sect. 2. Results of 
radio observations are presented in Sect. 3, while in Sect. 4 
we discuss radial profiles of the ICM surface brightness dis- 
tribution together with the geometry of the ICM. We de- 
scribe the search of the X-ray data for signatures of radio 
galaxy/ICM interactions in Sect. 5, and examine the pres- 
sure balance between the plasma of 3C 129 and the ICM in 
Sect. 6. Finally, we summarize our results in Sect. 7. 

We use H = h 100 km s" 1 Mpc" 1 with h = 0.65 
throughout, setting the cluster at a distance of 100 Mpc; 1 
arcsec corresponds to 474 pc. We quote all errors on 90% 
confidence level. All coordinates are given for equinox 2000. 
Distances from the radio core of 3C 129 are given as point 
to point distances, and not as distances along the radio tail. 

2 DATA SETS AND DATA ANALYSIS 

HI Observations with the Dominion Radio 
Astrophysical Observatory 

Since the 3C 129 area lies close to the Galactic plane it is in- 
cluded in the 21 cm line and continuum observations planned 
as part of the second phase of the Canadian Galactic Plane 
Survey (CGPS; see Taylor 1999 for an overview of its goals). 
The DRAO Synthesis telescope HI observations of this part 
of the sky were completed early in the second phase of the 
survey. They were reduced following the standard HI reduc- 
tion procedure for the CGPS, described by Higgs (1999). 
This work includes the insertion of "short spacing data" 
collected with the DRAO 26 m telescope and corrected for 
stray radiation effects. So the HI data are fully sampled out 
to an interferometer spacing of approximately 2900 wave- 
lengths. The resulting images have a resolution of 1 arcmin 
(Full Width Half Maximum, FWHM) in right ascension and 
1 arcmin cosec(<5) in declination 5. We used a method sim- 
ilar to that described by Normandeau (1999), to determine 
the optical depth in the direction toward the 3C 129 cluster. 

* The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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Briefly, HI spectra in the direction of the heads of the radio 
galaxies 3C 129 and 3C 129.1, which show absorption, were 
compared with spectra toward nearby "off source" positions. 
The comparisons indicate that the optical depth along the 
line of sight to these two strong sources averages to less than 
1. We assume that this optically thin condition holds any- 
where in the direction toward the 3C 129 cluster, so we can 
derive the HI column density from the expression: 



JVh = 1.823 x 10 18 x ZiT b (vi)Avi 



(1) 



Here, Nh is given in atoms cm -2 , the sum runs over all 
velocities Vi which contain Galactic hydrogen, Tb is the ra- 
dio brightness temperature in degree Kelvin, and Avi is the 
width of the ith velocity interval in km s -1 . 



Radio Observations with the Very Large Array 

The 330 MHz VLA data were acquired during 175 snapshot 
observations made over several years in the late 1990's. The 
data were taken in A, B, and C configurations. The final im- 
age has « 8"resolution (FWHM). The Root Mean Square 
(RMS) noise in the map is ~ 0.75 mjy beam - , and the 
Dynamic Range (DR, the signal-to-noise ratio of the peak 
brightness) is ~ 400. A detailed description of these obser- 
vations has been given by Lane et al. (2002). 

Data for the 1400 MHz image were obtained at the 
VLA on August 16, 2001 (4 hrs, C-configuration), and De- 
cember 6, 2001 (1.5 hrs, DnC-Configuration). Short obser- 
vations of 3C 48 were used for amplitude calibration, and 
the data were combined and mapped with the package MPS 
(Astronomical Imaging Processing System) using standard 
imaging techniques. The image has a FWHM resolution of 
« 18" x 14" at a position angle of -15°, an RMS noise of 100 
//Jy beam -1 , and a DR~ 3000. The primary aim of the VLA 
observations at L-band was to derive accurate 0.33 GHz - 1.4 
GHz spectral indices. The 330 MHz data have a minimum 
UV-spacing of 40A, and a maximum of 40.65kA (k stands 
for kilo); for the 1400 MHz data, the minimum spacing is 
180A and the maximum is 16.5kA. To make the spectral in- 
dex maps, both datasets were limited to a maximum UV 
of 15kA and a Gaussian taper was applied during imaging 
to match the measured beam sizes as closely as possible. A 
circular restoring beam of 18" x 18" was used in the final 
maps. We estimated the errors on the spectral indices by 
adding/subtracting 5 sigma to fluxes at one frequency and 
subtracting/adding 5 sigma from the other. This method 
gives conservative 5 sigma errors of 0.01 near the radio core 
and 0.25 at the end of the radio tails. 

Finally we used several VLA maps of 3C 129 and 3C 
129.1 with various effective angular resolutions that G. Tay- 
lor kindly provided. The pressure determination of the 3C 
129 radio plasma uses a 5 - 8 GHz spectral index map ex- 
tracted from two data sets with identical angular resolutions 
of 1.8"(Taylor et al. 2001). The width of the 3C 129 tail was 
determined at distances from 0"to 20", 0.3'to 2.5', and 2.5'to 
15' from the radio core from the 8 GHz, 0.83" angular resolu- 
tion map from Harris et al. (2002), the 8 GHz, 1.8" angular 
resolution map from Taylor et al. (2001), and the 0.33 GHz, 
8" angular resolution map from Lane et al. (2002), respec- 
tively. 

A high angular resolution map of 3C 129.1 (0.81" x 



0.76" FWHM) has been used for the radio contours shown 
in Fig. 5. 

Chandra Observations 

Our data consist of two pointings: a 31.4 ksec exposure taken 
on December 9, 2000 with the ACIS S, 6 chip array (aim- 
point close to 3C 129), and a 9.6 ksec exposure with the 
standard ACIS I CCD on January 9, 2001 (aimpoint close 
to 3C 129.1). The analysis uses the CIA0 2.2 software. Fil- 
tering for major background flares reduces the ACIS S data 
set to 28.3 ksec and the ACIS I data set to 9.3 ksec. We 
produce flat-fielded X-ray brightness maps using a soft (0.5- 
1 keV), an intermediate (1-2 keV) and a hard (2-5 keV) 
energy band. The brightness maps were corrected for in- 
strumental parameters such as the effective area etc. by di- 
viding by the appropriate exposure maps. The flat-fielded 
maps give the flux of 0.5-5 keV photons per cm 2 sarcsec 2 . 
We combine the flat-fielded maps of both observations by 
weighting each one with the corresponding integration time. 
Further details concerning these data can be found in Paper 
I. 

Some X-ray features have only modest statistical signifi- 
cance. Using the number of counts in the source "ON region" 
Aon, and the number of counts in the comparison "OFF re- 
gion" iVoFF, we compute the statistical significances based 
on the method of Li & Ma (1983). The method is based on 
Poisson probability distributions and correctly accounts for 
the statistical uncertainty of the counts in both the source 
and comparison regions. The statistical significance of an 
excess is given in standard deviations by: 

S*(iVoN,iVoFF;C) = V- 2 ln ( A ) ( 2 ) 
where 



A = 



f + 



and £ given by: 



Aon 
Won + Nqt, 

AqfF 



(3) 



f A(n)dQ 



A(Q) dn 



(4) 



where AfioN and AOoff are the solid angle areas of the 
source and background regions, respectively, and A(fl) is the 
the effective area of the detector. We use the exposure map 
produced for 1.7 keV photons as A(Q) and integrate over 
this map with the CIAO tool dmstat. Although the method 
is based on Poisson probability distributions, the expression 
for the statistical significance is only valid for Aon > 10 
and Aoff > 10 (Alexandreas et al. 1992). If this condition 
is not met, we use the equation of Alexandreas et al. (1992) 



2— J (1 + ?y) iV OFF+™+l 



(Aqff +n)\ 
JVqff! n\ 



(5) 



to compute the chance probability P c that Aon counts are 
observed by chance, taking into account the statistical un- 
certainty of the "OFF" measurement. We state the signifi- 
cances without taking into account "trial factors" due to the 
search for features over an extended solid angle region. 
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3 RADIO MAPS 

The N H -map derived from the HI observations, smoothed 
with a Gaussian of 2.45' FWHM is shown in Fig. 1. Although 
the N H -values change by 25% over the 1.25° x2. 5° field of 
view, the Chandra observations lie in a region of rather con- 
stant values of between 8.4xl0 21 cm" 2 and 9.2xl0 21 cm" 2 . 
Strasser and Taylor (2002) studied HI absorption spectra of 
273 background sources from the Canadian Galactic Plane 
Survey. In their Figure 2 they show that in optically thin re- 
gions the apparent off-source temperatures cover a range of 
about 100 K (upper bound) to 50 K (lower bound). Assum- 
ing a mean value of 75 K, these off-source data sets imply 
a maximum error in the HI determination of about 30%. 
Our results are marginally consistent with those of Leahy & 
Yin (2000) and Taylor et al. (2001), derived from fitting the 
cluster's X-ray energy spectra: the first authors report Nh = 
(5.7±0.3) x 10 21 cm" 2 based on ROSAT and EXOSAT data, 
and the latter authors get N H = (6.4±0.4) x 10 21 cm" 2 from 
ASCA data. 

The 1.4 GHz radio surface brightness distribution of the 
3C 129.1 and 3C 129 region is shown in Fig. 2, upper panel. 
The map shows substantially more detail than the previ- 
ously published map at this frequency (Jagers 1987). One 
clearly recognizes the peculiar radio feature (called "radio 
crosspiece" in the following) to the northeast of the radio 
head seen in earlier low-frequency 0.33 GHz data (Jagers 
& de Grijp 1983; Lane et al. 2002). Figure 2, lower panel, 
shows the 0.33 GHz-1.4 GHz spectral indices (a defined by 
S v oc v~ a ) derived from the 0.33 GHz data of Lane et al. 
(2002) and our 1.4 GHz data. The spectrum is flattest at the 
3C 129 radio core (a coro — 0.5). Between 0' and 6' from the 
radio core the spectrum of both radio tails stays rather con- 
stant with spectral indices from 0.5 to 0.8; between 6' and 
11' from the radio core it steepens to a w 1; and from 11' to 
15' it varies from 1 to 1.5. The region of the radio crosspiece 
has spectral indices of between 1.1 and 1.5. A natural inter- 
pretation of the spectral index distribution in the tail is that 
the electron population has an energy spectrum (dN/dE oc 
_E" P ) with p « 2.0. Within a distance of 8' from the radio 
core the radio spectrum steepens from a = 0.5 to 1.0, and 
thus the electron spectrum steepens from p — 2 to 3, as 
a consequence of synchrotron and Inverse Compton cooling 
of the radio plasma. Further cooling eventually results in 
a >1, once the electron spectra have cooled enough so that 
the high energy cutoff of electrons is low enough to affect the 
frequency band being observed. Weakening magnetic fields 
further down the radio tail will also contribute to this effect. 

Assuming no re-acceleration of electrons in the radio tail 
and using an order of magnitude magnetic field strength of 
B = 5/iG, we estimate the age of the electrons at 8' from 
the radio core where the spectral index of the synchrotron 
radiation is 0.5. For this purpose, we assume an initial elec- 
tron energy spectrum dN/dE oc _E" P with p = 2.0, and 
compute the 0.33 GHz-1.4 GHz spectral index taking into 
account synchrotron and Inverse Compton energy losses of 
the electrons (see Myers & Spangler, 1985). Following Jaffe 
& Perola (1974) we assume that electrons are isotropized on 
time scales short compared with their radiative lifetimes. In 
this approximation, Inverse Compton cooling on the Cos- 
mic Microwave Background can be take n into accoun t by 
using an effective magnetic field B' — B 2 + B^, MB with 



Bcmb = 3.25 (1 + z) 2 fiG to compute the energy losses of 
the electrons. The model reproduces the observed spectral 
index a — 1 for an age of 95 Myr of the electrons. This age 
translates into a velocity of the galaxy 3C 129 of: 

w g ai « 2260sin((9)- 1 (ho/0.65)" 1 km s" 1 (6) 

where 8 is the angle between the plane of motion of 3C 129 
and the line of sight toward the observer. Although the 
equipartition magnetic field at 8' from the radio core de- 
rived from the radio data is smaller (between 1.4 and 3.7 
/iG for lower frequency cutoffs v\ — 330 MHz and v\ — 
1 MHz, respectively) the radio plasma - ICM pressure bal- 
ance argues for a stronger magnetic field, possibly up to 15 
HG. Since v ga \ depends strongly on B, a range of allowed 
Wgai -values tightly constrains the magnetic field (see Sect. 
7). 

Lane et al. (2002) interpreted the radio crosspiece as 
fossil non-thermal electrons re-energized by the bow-shock 
of the radio galaxy 3C 129, supersonically traversing the 
ICM. Assuming that the segment of the energy spectrum of 
the electrons responsible for the 0.33 to 1.4 GHz emission in 
the radio crosspiece has not yet steepened due to radiative 
cooling, the spectral index a w 1.3 translates into a particle 
spectral index of p — 3.6. Invoking the standard theory of 
diffusive particle acceleration at shocks (Drury 1983), we 
derive an estimate of the Mach number of M ga i = 1.87. For 
an ICM temperature of 5 keV, the ICM sound speed is 1160 
km s" 1 , and M ga i = 1.87 corresponds to u ga i = 2170 km s" 1 
in fair agreement with above estimate. The difference of the 
radial velocities of 3C 129.1 [z = 0.022265) and 3C 129 (z = 
0.020814) is only 435 km s _1 , suggesting that the galaxy's 
velocity vector lies mainly in the plane of the sky. 

The spectral index values and variations found across 
the source between data at 0.33 and 1.4 GHz are in reason- 
able agreement with those found between 0.612 and 1.4 GHz 
(van Breugels 1982; Jagers 1987), and between 1.4 and 4.1 
GHz (Miley 1973). Comparing our results with the 0.074 
- 0.61 GHz spectral indices of Perley & Erickson (1979), 
we find that the spectrum flattens at low frequencies, while 
comparisons with the 0.408-2.7 GHz spectral indices of Riley 
(1973) and the 2.7-4.8 GHz and 4.8-10.6 GHz spectral in- 
dices of Feretti et al. (1998) show that the energy spectrum 
steepens at higher frequencies. 



4 THE GEOMETRY OF THE INTRACLUSTER 
MEDIUM AND DISCRETE SOURCES 

The upper panel of Fig. 3 shows the flat-fielded data of the 
two Chandra observations, summed with a weight propor- 
tional to the integration time and smoothed with a Gaussian 
of 18.8"FWHM. Neither before nor after smoothing the X- 
ray surface brightness do we detect obvious sharp gradients, 
like those found for a large number of more regular clus- 
ters (Markevitch et al. 2002) which might be produced by 
sloshing of the central ICM. The cluster has a considerable 
ellipticity. Using surface brightness iso-contours we derive el- 
lipticities e = (a — b)/a of about 0.3, a being the semi- major 
axis, and b the semi-minor axis. The surface brightness iso- 
contour ellipses are roughly aligned in the east-west direc- 
tion which depends slightly on the diameter of the ellipses. 
The surface brightness peaks at and around the location of 
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two discrete sources: an extended source near the core of 
the radio galaxy 3C 129.1 (a = 04 h 50 m 06 s .6308 and 5 = 
45°03'05.990"), and a point source, called FS1 in the follow- 
ing, at 2.0 arcmin toward the southwest of the center. We 
use the location a = 04 h 50 m 02 s .9 and 8 = 45° 02'26.0" that 
lies between these two sources as the cluster center. This lo- 
cation was determined by estimating centers of ellipses at 
various surface brightness levels. 

Another strong point source, called FS2 in the follow- 
ing, is found 9.2' from the cluster core toward its northwest. 
This second source had already been detected in ROSAT 
data by Leahy & Yin (2000). The location, X-ray flux, and 
information about radio counterparts of FS1 and FS2 are 
given in Table 3. Lacking optical spectra, we can not deter- 
mine whether these two source belong to the cluster or not. 
The only other known radio galaxy belonging to the cluster, 
WEIN 048 (Nilsson et al. 2000), was not in the field of view 
of the two Chandra observations. 

We searched for large scale irregularities of the ICM 
surface brightness by constructing models of the undis- 
turbed ICM surface brightness and subtracting these model 
from the surface brightness map. We used two models of 
the undisturbed ICM. The first model is a "standard King 
model" that is symmetric around the cluster center: 

=M - ^— — (7) 

[1 + (r/r c )f 2 

with the surface brightness S(r), the normalization constant 
So, the distance from the cluster center r, the cluster core 
radius r c , and the exponent /3. using the fiducial value /3 — 
2/3 we fit a core radius of r c = 7.1'with a \ -value of 317 
for 41 degrees of freedom. Fitting also the /3-parameter gives 
r c = 7.3' and ft = 0.68, and does not reduce the Rvalue 
substantially: 315 for 40 degrees of freedom. The surface 
brightness distributions after subtraction of this King model 
will show that the ICM distribution is quite irregular; we 
thus would not expect that a symmetric King model gives a 
good fit. 

The second model allows for variations of the ICM pro- 
file with polar angle. This model assimilates some of the 
large scale irregularities of the ICM and allows us to search 
for irregularities on smaller angular scales. Dividing the clus- 
ter into 18 pie sections of 20° opening angle each, 18 King 
models are fitted to the surface brightness. As shown in 
Fig. 4 for two example pie sections, the fiducial value of 
/3 = 2/3 yields unsatisfactory results: the observed surface 
brightness profiles show a more pronounced turn-down at 
cluster core distances between 5' and 10' than expected for 
[3 = 2/3. Leaving j3 as a free parameter, yields more sat- 
isfactory fits. For polar angle values (f> that are odd integer 
multiples of 10°, the surface brightness model is defined by 
the fitted King profiles. For intermediate values, we inter- 
polate between the predicted surface brightness of the adja- 
cent two King models. To avoid modification of the model 
by the radio galaxy 3C 129, we do not use the King model of 
the corresponding pie section (polar angle 240° -260°), but 
interpolate between the King models of the two adjacent 
pie sections. We use these surface brightness models only to 
find excess emission; the statistical significance of excess fea- 
tures on small angular scales is rigorously evaluated based 
on count statistics as described above. 



The residual surface brightness distributions smoothed 
with a Gaussian of 18.8"FWHM are shown in Fig. 3, cen- 
ter and lower panels. The map derived from subtracting the 
symmetric King model emphasizes the east-west elongation 
of the ICM: two regions with excess emission can be recog- 
nized: one is located at and toward the east of the cluster 
center, extending over at least 7'. The second one is located 
west/south- west of the cluster center. Note that the relative 
strength of these two excess regions depends strongly on the 
choice of the location at which we center the symmetric King 
profile. Moving this location by 1.3'toward the west of the 
cluster, reduces the western excess substantially, especially 
the bright region around the head of the radio galaxy 3C 129. 
Based on this surface brightness distribution, it is tempting 
to speculate that the excess region near the cluster center 
shows remnants of a sub-cluster passage through the cluster 
center. The curvature of the 3C 129 radio tail may be caused 
by a sub-cluster that attracts the radio galaxy and whose 
dark matter center of gravity even lies in the south-west of 
the south-western excess region. Upcoming observations of 
3C 129 with the XMM-Newton satellite will shed more light 
on the nature of these excess regions. The map obtained by 
subtracting the 18-slice pie-model from the surface bright- 
ness distribution does not show evidence for very sharp sur- 
face brightness gradients from large scale structure shocks or 
adjacent ICM phases of widely different X-ray emissivities. 



5 X-RAY SURFACE BRIGHTNESS NEAR THE 
TWO RADIO GALAXIES 

X-Ray Emission from 3C 129.1 

While the aimpoint of the ACIS I observation is located near 
3C 129.1, that of the ACIS S observation lies near the head 
of 3C 129 at about 12' from 3C 129.1. Thus, although the 
latter observation has a three times longer integration time, 
its angular resolution is substantially lower: 0.84"compared 
to 9"(FWHM) at the aimpoint of the observation and at 
12' off-axis, respectively. We therefore use only the ACIS I 
observation for the study of localized excess emission in the 
surrounding of 3C 129.1. 

The X-ray emission near 3C 129.1 is shown in Fig. 5 
together with a contour plot of the high resolution radio 
data. The high resolution radio map reveals a well defined 
northern jet, suggesting that the northern side approaches 
the observer. An excess of X-ray emission over that from the 
ICM can be recognized up to 10" from the radio core. The 
excess emission corresponds to a luminosity of 1.5xl0 41 erg 
s~ . There is clear evidence that at least some of this emis- 
sion comes from localized sources rather than from diffuse 
gas. At about 4" to the north of the radio core (source "C" ) 
we detect a total 7 counts (6 counts) from a sky region of 
1.1" (0.4") radius , where 0.25 counts (0.04 counts) had been 
expected based on an annular comparison region of ll"inner 
radius and 20"outer radius. In Table 4 we list the locations, 
statistical chance probabilities, and fluxes for the extended 
"center-source" as well as for 5 potential point sources. Using 
the Galactic neutral Hydrogen column density, we estimate 
a flux of 6.9xl0 40 erg s _1 for the center-source and fluxes 
between 8.7xl0 39 and 2.1 xlO 40 erg s _1 for the point source 
candidates. The 0.5-5 keV energy spectrum of the combined 
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excess X-ray photons can be described by a power law of 



spectral index a = 1.7±g;| (with EdN/dE cx E~ a ). Us- Bmin = (6 n) 2/7 (1+fc) 



ing only the events from the point source candidates, the 
best-fit spectral index is a — 1.3+° ®. Comparing the mean 
source spectrum luminosity with that of resolved sources in 
other X-ray faint galaxies (Sarazin, Irwin & Bregman 2000; 
Blanton et al. 2001a; Irwin, Sarazin & Bregman 2002), the 
luminosity of the 3C 129.1 sources seems to be relatively 
high, while their energy spectrum is relatively soft. This is 
especially true for source C which has the highest statistical 
significance, the highest flux and a very soft energy spectrum 
(see Table 4). 



X-Ray Surface Brightness Modification Near 3C 
129 

The evidence for X-ray emission from the 3C 129 radio core 
and jet has been described in Paper II. The excess emission 
was well localized and came from a sky region a couple of 
arc-seconds in diameter, located near the 3C 129 radio core. 
In Figs. 3 and 6, we excluded this core and jet emission 
from the analysis, by setting the residual surface brightness 
of these sky regions to zero. The residual X-ray map of Fig. 
3 (lower panel) and the enlarged view of the sky region near 
the radio head of 3C 129 in Fig. 6 reveal some evidence for 
an excess from a 12" x 36" (6 kpc x 17 kpc) large region at 
about 30" (14 kpc) toward the east of the radio core. Using 
an annulus of 45"inner radius and 90"outer radius centered 
on the radio core as comparison region, we estimate the 
statistical significance of the excess to be 3.7 a (ON: 131 
counts, OFF: 92 counts). 

The residual surface brightness map of Fig. 3 (lower 
panel) does not show any evidence for excess or deficit emis- 
sion along the tail of the radio galaxy 3C 129. Indeed, in ab- 
sence of strongly emitting sheets of compressed plasma, we 
do not expect to detect any excess or deficit, given the small 
volume of the tail. Using a rectangular test region with a 
solid angle coverage of 1.5' x 6' we can determine the ICM 
surface brightness to a la accuracy of between 4% and 15% 
near the radio head and near the extreme end of the radio 
tail, respectively. Assuming that (i) the tail is fully depleted 
of ICM, (ii) it lies in the plane of the sky, (iii) its depth along 
the line of sight equals its width perpendicular to the major 
tail extension, and (iv) that the ICM can be described by a 
King model with r c = 7.1'and (3 = 2/3 up to a cluster core 
distance of 40', we expect a reduction of the ICM surface 
brightness of only between 3% and 4% for locations near 
the radio head and near the extreme end of the radio tail, 
respectively. Although r c depends strongly on the location 
chosen as cluster center, the inferred reductions in the ICM 
surface brightness do not. 



6 PRESSURE BALANCE ALONG THE RADIO 
TAIL OF 3C 129 

In this section we compare the thermal pressure of the ISM 
and ICM with the minimum non-thermal pressure of the 
3C 129 radio tail. The energy density of the synchrotron 
plasma (including magnetic field energy) is minimized by 
the magnetic field (see Pacholczyk 1970): 




The fe-factor is the contribution to the non-thermal pres- 
sure from relativistic protons over and above that from the 
relativistic electrons, 4> is the volume filling factor for the 
radio emitting plasma, and V is the emitting volume, well 
constrained by the observations. The sum runs over the con- 
sidered frequency intervals i of the observed radio spectrum, 
the value of Ci2,i depends weakly on the synchrotron spec- 
tral index and the lower and higher cutoffs of the i th fre- 
quency interval, and L s ,i is the corresponding synchrotron 
luminosity. We compute true minimum energy estimates by 
assuming a magnetic field of B m i n , a volume filling factor 
of <f> = 1, no pressure contribution from electrons outside 
the energy range probed by the radio data, and no pressure 
contribution from relativistic protons [k =0). The finite 
angular resolution of the radio observations might overesti- 
mate the emission volume near the cores of the two radio 
galaxies and thus further lower the pressure estimate. The 
non-thermal particle energy density is given by: 



4 B 



3 8tt 



(9) 



and the equation for the minimum isotropic pressure of 
the synchrotron plasma (including magnetic field pressure) 
reads: 



, 2 -Bmin 

3 8tt 



(10) 



In the latter equation we conservatively assumed a tangled 
magnetic field. 

Based on these standard assumptions, Pmin can be com- 
puted from the radio data alone. As a starting point we 
use either the 0.33 GHz, the 5 GHz, or the 8 GHz surface 
brightness radio map. For each point along the radio tail 
we compute P m i n using thin cylindrical segments, in other 
words: disks, seen edge on. The disk-axes are aligned with 
the radio tail. The width (and depth) of a disk is given by 
the two points at which the radio intensity drops to 10% of 
its peak value. For each disk we compute the mean radio 
flux So at the frequency fo of the considered map, by av- 
eraging the surface brightness over the cross-section area of 
the disk. Given So, we construct the 0.33-8 GHz radio en- 
ergy spectrum, assuming the measured 0.33-1.4 GHz and 5-8 
GHz spectral indices over the frequency ranges from 0.33- 
2.65 GHz and 2.65-8 GHz, respectively. We use the ICM 
pressure from the spectroscopic ICM deprojection analysis 
of the Chandra data (Paper I) and neglect the pressure con- 
tribution of the ICM magnetic field. Taylor et al. (2001) 
estimate that the ICM magnetic field might be as strong as 
10 fiG. If this is true, the thermal and magnetic ICM pres- 
sures would be comparable at and near the head of the radio 
galaxy 3C 129. 

The pressure profiles are shown in Fig. 7. Over the first 
13" (6.1 kpc) from the nucleus of 3C 129 the minimum pres- 
sure of the radio plasma exceeds the pressure of the ICM 
by a factor of up to ~5, suggesting either self-confinement, 
or, confinement by the cold disk-like remnant of the ISM 
that has been stripped off the galaxy by the relative mo- 
tion of galaxy and cluster gas (Toniazzo & Schindler 2001). 
At 2' (57 kpc) from the radio core, the minimum pressure 
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fell by more than one order of magnitude and shows strong 
wiggles, suggesting "stormy local ICM weather". Unfortu- 
nately, low X-ray photon statistics do not allow us to deter- 
mine the ICM properties with sufficient statistical accuracy 
to assess the imprint of this stormy weather in the X-ray 
surface brightness or X-ray energy spectra. Further away 
than 2'from the radio core, the minimum pressure levels off 
at roughly l/40th of the ICM pressure. 

The equipartition magnetic fields lie between 16 pG and 
1 /LtG near the beginning and the end of the tail, respec- 
tively. The lower frequency cutoff of 330 MHz corresponds 
to Lorentz factors of 71 ~ 4000 (B/5 piG)~ 1/ ' 2 . We estimate 
the possible pressure contribution of electrons with very low 
Lorentz factors by extrapolating the radio energy spectrum 
down to v\ — 1 MHz, corresponding to Lorentz factors of 71 
~ 220 (B/5 fiG)^ 1 ^ 2 . Using the lower frequency cutoff, the 
inferred equipartition magnetic field strength increases to 
between 22 pG and 3 /iG near the beginning and the end of 
the tail, respectively. The corresponding pressure estimates 
are shown by the dashed dotted lines in Fig. 7. 

For the first 10' of the tail, the radio spectrum is suffi- 
ciently flat so that the lower cutoff frequency does not quali- 
tatively alter the result: the minimum non-thermal pressure 
lies substantially below the ICM pressure. We will discuss 
the implications in the next section. 



7 SUMMARY AND DISCUSSION 

The Chandra observations show that the radio galaxy 3C 
129.1 is associated with an extended X-ray excess. The ex- 
cess can be divided into 5 point sources and a central com- 
ponent. On the small angular scale of the central compo- 
nent (4"diameter) the ISM is not expected to contribute to 
the excess. We therefore consider it most probable that the 
central excess emission is produced by galactic sources and 
possibly also by the central AGN. 

For the radio galaxy 3C 129 we find, in addition to 
the core and jet emission described in Paper II, some evi- 
dence for extended X-ray emission from a 12" x 36" (6 kpc 
x 17 kpc) large region to the east of the radio core. If con- 
firmed with higher statistical significance, a natural expla- 
nation of the excess would be thermal emission from plasma 
compressed and heated by the bow shock that precedes the 
galaxy. 

Comparing the minimum pressure of 3C 129 radio tail 
with the thermal ICM pressure, the radio plasma appears 
to be under-pressured by a factor of between 10 and 40. 
Over the first 10'of the radio tail, the radio spectrum is 
sufficiently flat so that the results depend only weakly on 
the assumed minimum Lorentz factor of the non-thermal 
electron distribution. 

Remarkably, a very similar pressure mismatch has been 
found for the kilo-parsec radio features of a number of FR I 
and FR II radio galaxies (see Hardcastle & Worral 2000 and 
references therein). For most radio galaxies the magnetic 
field inside the radio plasma is poorly constrained toward 
higher values and might dominate the plasma pressure. In 
the case of 3C 129, we have an additional handle on the 
magnetic field strength from the constraints on the radia- 
tive age of the radio-electrons: If the missing pressure of the 
radio plasma is provided by the magnetic field alone, the 



implied field strengths of between 15/xG and 30pG trans- 
late into very high lower limits on the galaxy velocity of be- 
tween 11,000 km s _1 and 31,000 km s _1 (sec Eq. 6). Buoyant 
motion of the radio plasma would proceed with a velocity 
lower or equal to the ICM sound speed, and would thus not 
qualitatively alter our conclusions. The apparent pressure 
mismatch thus suggests one of the following possibilities: (i) 
thermal or relativistic protons and/or thermal or relativis- 
tic low-energy electrons carry the major fraction of the radio 
plasma's pressure; (ii) the volume filling factor of the radio 
plasma <f> is substantially smaller than unity: (f> = 0.0016 - 
0.018; (iii) the B-field is much larger than the equipartition 
magnetic field and electrons are re-accelerated inside the ra- 
dio tail; (iv) the radio galaxy lies behind or in front of the 
plane of the cluster so that the true ambient ICM pressure 
is lower than assumed in our analysis. The detection of an 
"ICM cavity" associated with the radio tail would make it 
possible to estimate the volume filling factor. We hope to 
detect such a cavity with an upcoming XMM-Newton ob- 
servation of 3C 129. 
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Table 1. 3C 129 Chandra Data Sets. 



Description 


Energy 


Integration Time 


Ang. Res. 


Reference 








(FWHM) 




Chandra ACIS I 


0.5-8 keV 


9.3 ksec 


0.84"(aimpoint) 


Krawczynski 2002 


Chandra ACIS S 


0.5-8 keV 


28.3 ksec 


0.84"(aimpoint) 


Krawczynski 2002 



Table 2. Radio Data Sets Used in this Paper. 



Frequency 


Instrument 


New Data 


Ang. Res. 
(FWHM) 


Noise RMS 
[mjybeam -1 


DR 

1 


Reference 


HI 


DRAO 


Y 


V 






this work 


330 MHz 


VLA 


N 


8" 


0.75 


400 


Lane et al. 2002 


1.4 GHz 


VLA 


Y 


16" 


0.1 


3000 


this work 


4.7 GHz (3C 129) 


VLA 


N 


1.8" 


0.072 


545 


Taylor et al. 2002 


8 GHz (3C 129) 


VLA 


N 


1.8" 


0.04 


925 


Taylor et al. 2002 


8 GHz (3C 129) 


VLA 


N 


0.83" 


0.03 


1287 


Harris et al. 2002 


8 GHz (3C 129.1) 


VLA 


N 


0.8" 


0.012 


225 


Taylor et al. 2002 



Table 3. X-ray sources not associated with known cluster galaxies. 



Name 


RA (Eq. 2000) 


Decl (Eq. 2000) 


0.5-5 keV Flux° 


Radio Counterpart 






[hh:mm:ss.ss] 


[dd:mm:ss.ss] 


[ergcm~ 2 s _1 j 






FS1 


04:49:58.5 


45:01:44.8 


5.13xl0~ 14 


0.24 mjy @4.8 GHz, 0.11 mjy @i 


3 GHz 


FS2 


04:49:54.8 


45:11:31.4 


8.85X10- 14 


no 





a Observed flux. 



Table 4. X-ray Sources of 3C 129.1 



Name 


RA (Eq. 2000) 


Decl (Eq. 2000) 




n N a 
'1 JV OFF 


Chance 




£ 0.5-5keV 


Flux c 


Luminosity d 




[hh:mm:ss.ss] 








Probab. 




[keV] 


r -2 -ii 

erg cm s 


[ergs-] 


A 


04:50:07.075 


45:03:02.91 


3 


0.40 


8.0x10" 


3 


2.17±0.29 


2.0 10~ 15 


8.1 xlO 39 


B 


04:50:06.729 


45:03:02.92 


4 


0.40 


7.9x10- 


-4 


3.58±0.68 


2.7 10" 16 


1.1 xlO 40 


C 


04:50:06.618 


45:03:09.18 


7 


0.40 


2.3x10" 


-7 


2.27±0.17 


5.0 10" 1S 


2.1 xlO 40 


Center 


04:50:06.584 


45:03:05.92 


25 


2.80 


6.7x10- 


-16 


2.48±0.15 


1.7 10- 14 


7.0X10 40 


D 


04:50:06.149 


45:03:01.35 


5 


0.40 


6.2x10- 


■5 


3.69±0.58 


3.5 10" 1S 


1.5X10 40 


E 


04:50:05.740 


45:03:03.61 


6 


0.40 


4.1x10- 


6 


4.16±0.46 


4.2 10" 16 


1.8 xlO 40 



a Expected number of "background" counts in the ON region. 
6 Mean energy of excess 0.5-5 keV Chandra photons. 
c Observed flux. 

d "Unabsorbed luminosity", assumes N H = 8.8xl0 21 cm -2 . 
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Galactic longitude 
161.200 160.800 160.400 160.000 




Figure 1. Neutral Hydrogen Column density in the direction of the cluster 3C 129 derived from the DRAO HI data, smoothed with a 
Gaussian of 2.45' FWHM. Darker regions show higher Njj column densities. Rather atypically, the column densities arc a bit lower along 
the galactic plane than at slightly lower and higher latitudes. The ten contour levels are linearly spaced from 8xl0 21 to 10 22 cm~ 2 . 
Compared to nearby fields in the Galactic plane (see e.g. Leahy and Roger 1991, Figs. 7 and 8), only modest gradients of the Ntr-column 
density arc observed. 
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Figure 2. The upper panel shows the 1.4 GHz surface brightness distribution. The peak surface brightness is 0.29 Jy beam -1 and the 
8 contours are logarithmically spaced from 0.001 to 0.256 Jy beam" 1 (bcamwidth 16"FWHM). Darker regions correspond to a higher 
surface brightness. The lower panel shows the 0.33-1.4 GHz spectral indices (a from S v oc v~ a ) derived from the 0.33 GHz data of Lane 
et al. (2002) and the 1.4 GHz data of the upper panel. Darker regions show relatively steeper energy spectra (larger a) and the 7 contour 
levels are linearly spaced from a = 0.5 to a = 2. (Aa = 0.25). The small radio source in the east is 3C 129.1 and the larger one in the 
west is 3C 129. 
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Figure 3. The gray-scale plot in the upper panel shows the flat-fielded Chandra 0.5-5 keV data (ACIS I and ACIS S), smoothed with 
a Gaussian of 18.8"FWHM. Segments of ACIS-I chip gaps are visible as white stripes outside the regions covered by the ACIS-S chips. 
The center panel shows the residual surface brightness distribution after subtraction of the symmetric ICM model; the lower panel shows 
the residual surface brightness after subtraction of the 18-slice pie model. In all three panels the gray-scales are linear and high emission 
regions are darker than low emission regions; the surface brightness in the upper panel peaks at the location of FS 1 with 3.8 xl0~ 6 
photons cm~ 2 s _1 arcsec~ 2 , at the radio core of 3C 129 it is l.OxlO -6 photons cm~ 2 s~ 1 arcsec -2 ; the gray of the "noise" at the very 
end of the 3C 129 tail is 10 — 7 photons cm~ 2 s~ 1 arcscc~ 2 ; The center and lower panels use the same gray-scales with values of 2.1xl0 — 6 
photons cm~ 2 s _1 arcscc~ 2 at the location of FS 1, and 10 — 7 photons cm~ 2 s~ 1 arcsec -2 at the noisy gray areas near the end of the 
tail of 3C 129, respectively. The radio contours have been derived from the 330 MHz radio map and indicate the location of the radio 
galaxy 3C 129.1 at the cluster center and 3C 129 to the east of the cluster center. Contours are drawn at 0.0025, 0.035, 0.08, 0.12, 0.15, 
and 0.20 Jy beam" 1 (beamwidth 8"FWHM). 
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Figure 4. Radial 0.5-5 keV surface brightness distributions. The left and right profiles have been derived from "pie-sections" of the 
Chandra data with polar angles from 80° to 100° and from 240° to 260° , respectively. The full lines show fits of King models and the 
best-fit values of the normalization factor So (in 10 -8 photons cm™ 2 s™ 1 arcsec -2 ), core radius r c (in arcmin), and /3 are given in the 
upper right corner of each panel as PI, P2, and P3, respectively. In the left hand panel, the dotted line shows the fitted profile of the 
right hand panel, and emphasizes the difference of the two profiles. In the right hand panel, the dotted line shows a King model with 
= 2/3. 
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Figure 5. The gray-scale plot shows the 0.5-5 keV X-ray surface brightness at the location of the radio galaxy 3C 129.1, smoothed 
with a Gaussian of 1.1"FWHM. The gray-scale is linear and peaks near the radio core of 3C 129.1 with 2.5xl0 — 7 photons cm -2 s —1 
arcsec™ 2 ; for source A the surface brightness peaks at lxlO™ 7 photons cm™ 2 s™ 1 arcsec -2 . The radio contours have been derived from 
the 8 GHz map with a beamwidth of 0.81" X 0.76" FWHM with 9 logarithmically spaced levels from 2 to 512 mjy beam™ 1 . The high 
resolution radio map reveals a well defined northern jet, suggesting that the northern side approaches the observer. The photon excesses 
discussed in the text have been marked by circles. 
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Figure 6. The gray-scale plot shows the residual 0.5-5 keV surface brightness distribution near the head of the radio galaxy 3C 129 after 
subtraction of the 18-slice pie model, smoothed with a 9.4"FWHM Gaussian. The gray-scale is linear and peaks near the radio core of 
3C 129 with a value of 7.4xl0 -7 photons cm 2 s 1 arcscc 2 ; the two bins at the bottom left side is a similar gray- tone have a surface 
brightness of 5.7xl0 -7 photons cm -2 s -1 arcscc -2 . The radio contours have been derived from the 5 GHz VLA data with 1.8"FWHM 
beamwidth and the contours are logarithmically spaced from 1 to 32 mjy beam -1 . An X-ray bright region southeast of the radio core 
can be recognized. The rectangle shows the region of the sky that we used to determine the statistical significance of the excess. 
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Figure 7. The thermal ICM and non-thermal radio plasma pressures for the radio galaxy 3C 129 as function of distance to the cluster 
center at a = 04 h 50 m 02 s .9 and <5 = 45° 02'26.0" (J2000). The thermal ICM pressure from the spectroscopic deprojection analysis of 
the X-ray data (Paper I) is shown by the solid line and the data points. The dashed and dotted lines show the minimum pressure of 
the radio plasma for the northern and southern radio tail, respectively, computed for a minimum radio frequency v\ = 330 MHz. An 
estimate of the pressure contribution from very low energy electrons is shown by the dashed dotted lines, which have been computed 
under the assumption that the radio spectrum extends down to v\ =1 MHz with the same spectral index as measured for 330 MHz - 
1.4 GHz. 
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